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A computer model was developed for estimation of the kinetics of microbial inactivation by pulsed 
electric field. The model is based on the electroporation theory of individual membrane damage, 
where spherical cell geometry and distribution of cell sizes are assumed. The variation of microbial 
cell sizes was assumed to follow a statistical probability distribution of the Gaussian type. Surviving 
kinetics was approximated by Weibull equation. The dependencies of two Weibull parameters (shape 
n and time t, respectively) versus electric field intensity E and width of cell diameters distribution 
was studied. 
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I. INTRODUCTION 

Pulsed electric fields (PEF) processing is a promising 
method of food preservation. Many investigators have 
shown the effectiveness of PEF application for killing 
bacteria in liquid foods (Barbosa-Canovas et al., 1998; 
Barsotti & Cheftel, 1998; Wouters & Smelt, 1997). How- 
ever, there still exist a considerable gap in understanding 
the inactivation mechanism. The important problem is 
to elucidate how kinetics of killing depends on the type of 
bacteria and treatment protocol (electric field strength E, 
form of pulse, pulse duration t^, total time of treatment). 

The phenomenon of PEF-inactivation of microbial 
cells is related to selective damage of biological mem- 
brane. Electrical conductivity of the membrane a is 
very low. The reported values are of order of a ~ 
lO^-lCT^^TO-^Kotnik et al., 1998). Therefore, the 
highest drop of potential occurs on the membranes. The 
transmembrane potential of a spherical cell u m depends 
on the angle between the external field E direction and 




FIG. 1: Spherical microbial cells in external field E. Cell di- 
ameters d c are assumed to follow a distribution function of 
the Gaussian type F(d c ). Here 6 is an angle between elec- 
tric field direction E and a radius vector r at the surface of 
membrane. 



the radius-vector r on the membrane surface, where po- 
tential is to be determined (Fig. 1). This potential may 
be determined using the well-known Schwan's equation 
(Schwan, 1957), 

u m = 0.7hfd c Ecos6, (1) 

where d c is the cell diameter, and / is a parameter de- 
pending on electrophysical and dimensional properties of 
the membrane, cell and surrounding media. In dilute 
suspension of cells, parameter / is close to 1 (Kotnik et 
al, 1998). 

So, the cell transmembrane potentials u m in suspen- 
sion of cells depend on the cell diameter d c and angle 
and for the single cell value of u m is maximal at cell 
poles and decreases to zero at 6 = ±tt/2. That is why the 
membrane damage probability is maximal at membrane 
poles and that bigger sized microbial cells are killed be- 
fore smaller ones. Microbial cells always show a variety 
of shapes and dimensions (Bergey, 1986). Their size may 
vary depending on their age, the nutrients in the growth 
medium, release mechanisms of microbial particles etc. 
(Harding, 1975; Reponen et al., 1992). So, microbial 
cells killing probability can also change from cell to cell. 

After PEF treatment during time t in the electric field 
E the surviving fraction S(t, E) is defined as the ratio 
of the number of undamaged microbial cells to the to- 
tal number of microbial cells (Barbosa-Canovas et al., 
1998). If all the cells are spherical and are of the same 
size, then their damage may be considered as statistically 
independent events, and the time dependence S(t, E) can 
be approximated by the first-order kinetic equation: 

S(t,E) = exp(-t/r(E)), (2) 

where t(E) is a time parameter that corresponds to the 
effective inactivation time of cells in the external electric 
field E. 

Unfortunately, in most cases, the simplest approxi- 
mation of the first-order kinetics is not applicable for 
description of the microbial inactivation experiments in 



2 



pulsed electric fields (Barbosa-Canovas et al., 1998). Hul- 
sheger et al. (1983) proposed an empirical equation of 
type 



S(t,E) = (t/t c ) 



-(E-E c )/k 



(3) 



where t c and E c are the threshold treatment time and 
electric field intensity, and k is an empirical parameter. 
Although this equation is very popular, it has no theo- 
retical justification. 

The models widely used presently for describing the 
survival curves are Fermi, log-log and log-logistic models 
(Barbosa-Canovas et al., 1998; Peleg, 1996; Alvarez et al., 
2000) but they are also of an empirical nature. Specifi- 
cally, the Weibull distribution may be a useful general- 
ization that that includes exponential first order kinetics 
as a special case (Peleg, 1995, 1999) 



S(t,E)=exp(-(t/r(E)r^), 



(4) 



where t{E) is a time parameter and n(E) is a shape 
parameter. In the case when n(E)=l, Eq. (4) reduces to 
Eq. (2). 

The time parameter t(E) in Weibull distribution ac- 
counts for the effective inactivation time, and shape pa- 
rameter n(E) accounts for the concavity of a survival 
curve (van Boekel, 2002). Weibull distribution was ap- 
plied for fitting experimental PEF inactivation data (Al- 
varez et al., 2002), but physical meaning of the obtained 
parameters t{E) and n(E) was not elucidated yet. 

A possible deviation from the first order PEF inacti- 
vation kinetics may be caused by existence of a variety 
of microbial shapes and dimensions. The purpose of this 
paper is to analyse how the form of the survivor curves 
can reflect existence of a distribution of cell diameters. 



II. COMPUTATIONAL MODEL AND DETAILS 
OF CALCULATIONS 

The Monte Carlo technique was used for simulation of 
PEF inactivation kinetics of microbial cells. Initial num- 
ber of microbial cell in suspension was put as N o =10 7 . A 
Gaussian law distribution function of cell diameters was 
assumed (Fig.l) 



F(d c ) 



2ttA 



exp(- 



(dc - d c f 
2A 2 



(5) 



where d c and A represent the average diameter and the 
standard deviation, respectively. 

An arbitrary microbial cell was chosen in suspension 
for a given time t. Then, a random point on the mem- 
brane surface was chosen by generating of cos# value ran- 
domly in the intervals — 1 — 1—1- The lifetime of a mem- 
brane r on the surface of a cell depends on its diameter 
d c , angle 9, and intensity of external field E. It was found 
on the basis of the transient aqueous pore model (Weaver 
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FIG. 2: Calculated survivor curves S(t/Too)- Insert shows the 
shape n and relative time t /too parameters of Weibull distri- 
bution versus the relative width of cell diameter distribution 
A/d c . All the calculations were made for E* = 10. 



& Chizmadzhev, 1996), that: 



r(9,d c ,E) = TooCXp 



iru) 2 /kTj 



l + (u m (9,d c ,E)/u )^ 



(6) 



where u m (d c , 9) was calculated from Eq. (1). Here, is 
the parameter (t — > in the limit of very high electric 
fields), to and 7 are the line and surface tensions of mem- 
brane, respectively, k is the Boltz mann constant, T is 
the absolute temperature, u = y/2j/(C m (£w/£m — 1)) 
is the voltage parameter (the dimension of u Q is Volts), 
C m is the specific capacitance of a membrane, s w , e m are 
the relative dielectric permittivities of the aqueous phase 
and of the membrane, respectively. 

The probability of the chosen cell damage was 
approximated by the first-order kinetic equation as 
exp(—t/T(9,d c ,E)). This procedure was repeated for all 
the cells in the suspension. Then, the number of killed 
cells was enumerated for the given time t, surviving frac- 
tion S(t, E) was calculated, time was increased by a time 
step and procedure was repeated from beginning. 

In this work, the voltage scale parameter was estimated 
as u a w 0.17V from data obtained by Lebedeva (1987) 
for the general lipid membranes (u) s=s 1.69 * 10~ n N, 
7 w 2*10~ 3 N/m, e w w 80, e m w 2, C m w 3.5*10~ 3 F/m 2 
at T — 298K). The time scale parameter was put as 
~ 3.7 * 10~ 7 s (Lebedeva, 1987). Dimensionlcss re- 
duced field intensity was defined as E* = E/E a , where 
E = u / (0.75 fd c ) was estimated as E a w 2.27 kV/cm 
at d c i=a 1/im. 



III. RESULTS AND DISCUSSION 

Figure 2 presents some examples of the calculated sur- 
vivor curves S(t) for suspension of cells (symbols) at the 
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FIG. 3: Shape n and relative time t/too parameters of 
Weibull distribution versus reduced electric field intensity 
E* = E(0.75fd c /u o ) for suspension of cells with diameter dis- 
tribution for different values of relative width of cells diameter 
distribution A/d c . The simulated data were fitted with the 
Weibull equation within time interval of < t/roc < 200. 



given value of reduced field intensity E* — 10 (E «22.7 
kV/cm, at d c = 1/im). The first order kinetics law is only 
observed for suspensions of identical cells (A/d c = 0). In 
other cases, the noticeable deviations from the first or- 
der kinetics are observed and the more pronounced devi- 
ations are observed with increase of A/d c . 

The solid lines drawn through the symbols are the best 
fit to data simulated using the Weibull function (Eq. 4). 
The Weibull law seems to be very appropriate for approx- 
imation of calculated survival curves (in all the cases the 
correlation coefficient p was higher than 0.993). Insert 
to Fig. 2 shows the shape n and relative time t/too pa- 
rameters of Weibull distribution versus relative width of 
cell diameter distribution A/d c . In these estimations the 
fitting was done within the time interval 1 < i/r^ < 200. 
At Too ~ 3.7 * 10~ 7 s (Lebedeva, 1987) this time interval 
corresponds to < t < 7Aps. Both the shape parame- 
ter n and the relative time t/too parameter initially de- 
crease with increases of the standard deviation A. Then, 
a small elevation of these values is observed, which can be 
explained by the distortion of the Gaussian distribution 
at higher values of A/d c . But this model always gives 
only upward concavity, i.e. n < 1. 

Parameters n and t/too are also very sensitive to the 
value of electric field intensity E (Fig. 3). Dependencies 
of n versus E* are rather complex, but in all cases param- 
eter n increases with increase of the distribution width 
A/d c . Relative time parameter t/t^ decreases consid- 
erably with field intensity E increase, but it is practi- 
cally independent of A/d c . So, it is possible to conclude 
that upward concavity of survivor curve is rather sensi- 
tive both to the field intensity and to the variability of 
microbial diameter distribution, but the effective inacti- 




FIG. 4: Shape parameter n and relative time parameter t/too 
versus time interval t ma i/7"oo. Calculations were done at re- 
duced electric field intensity of E* — E(0.75fd c /u o ) = 6 and 
relative width of cell diameter distribution A/d c — 0.2. Dash 
line shows the time t m ax /t^ — 200 used for calculation of 
data presented in Fig.l and Fig. 2. The insert shows cal- 
culated survivor curve S(t/Too) for this particular case. The 
solid line is the best fit to the simulated data (symbols) with 
the Weibull equation in the time interval i/roo < 10 5 . 



vation time t(E) is insensitive to the variations in cell 
diameters. 

The numerically estimated shape n and relative time 
t/tqo parameters are rather sensitive to the time interval 
of Monte Carlo data fitting with Weibull function Eq. 
(4). Fig. 4 shows a typical example of n and t/t^ versus 
tmax/Too dependencies for fitting of the same survival 
curve when upper bound of the time interval 1 < t/roo < 
tmax/Too is varied. In fact, t max /Too is a relative total 
time of treatment. In all cases, the apparent consistency 
between Monte Carlo data and Weibull function with 
adjusted parameters n and t/too was rather good, and 
the correlation coefficients p lied in the interval 0.993- 
0.998. But at the same time, parameter n and t/too are 
sensitive to the upper cutting boundary t max /Too, and 
this fact reflects existence of an intrinsic inconsistency 
between unknown survival function and Weibull function. 



IV. CONCLUSION 

The discussed illustrative examples show that geome- 
try of the survival curve is very sensitive to the distri- 
bution of cell diameters. The Weibull function seems to 
be appropriate for approximation of the calculated sur- 
vival curves. The parameters of this function r and n 
are rather sensitive to the width of distribution of the 
cells diameters, electric field intensity and total time of 
treatment. We would note, that the proposed model is 
based on several restrictive assumptions. The survival 
kinetics may reflect many intrinsic details of the real mi- 
crobial cells. For explanation of the survival curves, ex- 
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pcrimentally observed for PEF-inactivated population, it 
is necessary to introduce into the model the experimen- 
tally determined distribution functions of cell diameters. 
It is also desirable to use in calculations more realistic 
law of an individual membrane damage, based on exper- 
imentally estimated data for the given bacterial popu- 
lation. The possible effects of sub-lethal damage, when 
bacterial damage needs some critical destructive expo- 
sure, also were not considered in this model. A deviation 
from the first-order kinetics may be also influenced by 
others factors, such as existence of bacterial geometry 
anisotropy and distribution of bacterial orientations. So, 
in future it is seems to be important to find correlations 
between variations in factors, influencing bacterial ge- 
ometry, dimension distribution function, details of mem- 
brane damage and parameters of bacterial inactivation 
kinetics. Such work should be done in order to improve 
practically important PEF-treatment regimes for reach- 
ing a desirable value of microbial inactivation. 
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